Impurity states and interlayer tunneling in high temperature superconductors 
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We argue that the Scanning Tunneling Microscope (STM) images of resonant states generated 
by doping Zn or Ni impurities into Cu-O planes of BSCCO are the result of quantum interference 
of the impurity signal coming from several distinct paths. The impurity image seen on the surface 
is greatly affected by interlayer tunneling matrix elements. We find that the optimal tunneling 
path between the STM tip and the metal (Cu, Zn, or Ni) d x 2_ y 2 orbitals in the Cu-O plane involves 
intermediate excited states. This tunneling path leads to the four-fold nonlocal filter of the impurity 
state in Cu-O plane that explains the experimental impurity spectra. Applications of the tunneling 
filter to the Cu vacancy defects and "direct" tunneling into Cu-O planes are also discussed. 
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Recently J.C. Davis and collaborators applied the STM 
technique to image single Zn and Ni impurities in opti- 
mally doped BSCCO These experiments proved 
that one can image single impurity states in an uncon- 
ventional superconductor and demonstrated the highly 
anisotropic structure of these states. 

Although it appears that on a gross scale these find- 
ings can be understood in terms of a conventional d-wave 
superconductor perturbed by potential scattering, upon 
closer inspection problems of principle seem to arise. The 
impurity states observed by STM are characterized by 
two main features: (1) energy and width of the impurity- 
induced resonance in the density of states (DOS), and (2) 
the spatial structure of the resonance. While the DOS 
seems to be satisfactorily described by a single-site im- 
purity model H, the real space distribution of intensity 
cannot be fit by this model. The main problem with the 
Zn impurity image seen in the STM experiments is that 
the intensity of the signal on the impurity site is very 
bright, which is at odds with the unitary scattering off 
Zn. We remind that Zn 2+ has a closed d shell and hence, 
exactly on the impurity site the scattering potential is 
very strong. Unitary scattering is equivalent to the hard 
wall condition for the conduction states and therefore no 
or very little intensity of electron states is expected on the 
Zn site. A similar problem arises also with explaining the 
Ni-induced resonance. 

Here we demonstrate that these problems find a nat- 
ural resolution in terms of the specific way in which the 
local density of states of the cuprate planes is probed 
in the STM experiments. We argue that the quantum- 
mechanical nature of the tunneling from the STM tip 
into the Cu-O layer that hosts impurity requires tunnel- 
ing through the uppermost insulating Bi-0 layer which 
effectively filters the signal. Surprisingly, similar filter- 
ing should also take place even in the case of "direct" 
tunneling into Cu-O plane. Such non-local tunneling has 
profound consequences for the real space image of the 
impurity state seen by STM. 

There are two major types of the tunneling routes be- 



tween the STM tip and the conducting orbitals in the 
Cu-O plane: (a) Direct tunneling due to the overlap be- 
tween the tip and the planar 3d x 2_ y 2 wavcf unctions, and 
(b) indirect tunneling through intermediate excited (oc- 
cupied or empty) states, Fig. |I|a. The direct tunneling 
probability over the experimentally relevant distances of 
about 10A is however exponentially small [Q. On 
the other hand, the importance of the excited states for 
mediating STM tunneling follows directly from the Bi-0 
topographs that clearly show the positions of the nomi- 
nally insulating Bi atoms on the surface of BSCCO. The 
analysis of the topographs implies that the excited Bi 
orbitals focus the flow of the tunneling electrons. There- 
fore, we argue that the indirect tunneling via overlapping 
intermediate orbitals is the dominant tunneling mecha- 
nism. 

The strongest indirect tunneling channel involves in- 
termediate states that have the largest overlaps. In the 
case of BSCCO, these are the orbitals that extend out 
of the planes, such as 4s or Zd^ z 2_ r 2 of Cu and 6p z of 
Bi. Furthermore, only the states with zero in-plane or- 
bital momentum have non-zero overlap with the apical 
oxygen 2p z and 3s orbitals that play important role in 
the interlayer communication. Being radially symmetric 
in the Cu-O plane, such states cannot couple to the rel- 
evant 3d x 2_ y 2 states on-site. They do couple however to 
the neighboring 3d x 2_ y 2 orbitals through the d- wave-like 
fork (Fig. ^]b). The resulting tunneling amplitude is 
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where vE^j is the impurity state wavefunction on site 
Hence, the symmetry analysis leads to the sur- 
prising conclusion: Tunneling on top of a particular Cu 
(or Zn or Ni) atom in the Cu-O plane does not probe 
its 3d x 2_ y 2 orbitals, but rather measures a linear com- 
bination of its neighbors. The tunneling matrix ele- 
ment determines the intensity of the impurity signal 
■^[^ij] = l-^j-jl 2 - Calculating the spectral intensities in 
the vicinity of the impurity using this filtering function 
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reproduces the features seen in the experiments in great 
detail (Figs. || and ||). The aforementioned anomalies 
find an explanation in terms of interferences associated 
with the non-local way in which the electronic states are 
probed. The described filtering mechanism is related to 
the one responsible for the interlayer tunneling in the 
bulk cuprates |^,^]. An alternative non-local filter based 
on the incoherent direct tunneling from the tip into Cu-0 
plane was recently studied by Zhu et al. (7| . The angular 
effect of the filter on the far asymptotic of the impurity 
states was discussed earlier S . 
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FIG. 1. (a) Tunneling from STM tip into 3d x 2_ y 2 orbitals 
in the Cu-O layer occurs through the virtual energy states 
of Bi and Cu. The tunneling processes involve both empty 
(top) and occupied (bottom) intermediate states. All such 
processes contribute coherently to the tunneling amplitude 
from the tip to the 3d x 2_ y 2 orbitals. (b) Tunneling through a 
particular Bi atom does not probe the 3d x 2_ y 2 orbital of the 
metal atom (Cu, Zn or Ni) right underneath it due to the van- 
ishing overlap. Instead, the tunneling involves orbitals that 
extend out of the planes, such as 4s (shown) or 3d 3z 2_ r 2 of 
Cu and 6p z of Bi. These orbitals are symmetric in the Cu-O 
plane and hence couple to the neighboring_met&\ 3d x 2_ y 2 or- 
bitals through the d- wave-like fork, Eq. (jlh . 



Similar to the Zn case, the Ni impurity tunneling inten- 
sity is maximal at the Bi position immediately above the 
impurity site. However, there are several important dif- 



ferences between the observed Zn- and Ni-induced states: 
(1) The energy of the Zn state is close to the chemical po- 
tential, e Zn = — 2 mV, while the Ni state energy is larger 
and is split, e Ni = 9 and 18 mV, (2) the Zn state appears 
only on the negative bias, while the Ni state shows up 
both on positive bias and the symmetric negative bias. 
In this letter we demonstrate that these experimental fea- 
tures can also be reproduced within the standard theory 
of the impurity states || , with the spatial structure of the 
states being reproduced by properly taking into account 
the fork effect. 

The starting point of our model is a two-dimensional 
mean-held (MF) Hamiltonian with the nearest neighbor 
attraction, V, which yields d-wave superconductivity in 
the range of dopings close to half-hlling, 
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Here, Ay = V(cnCj^) is the self-consistent MF super- 
conducting order parameter. The hopping tij equals t 
for nearest neighbors and t' for for the second-nearest 
neighbor sites i and j. The parameters that are relevant 
for BSCCO are t = 400 meV, t' = -0.3t. To match the 
amplitude of the superconducting gap in optimally-doped 
BSCCO, which is about 40 meV, we choose the attrac- 
tion V = — 0.525i. The chemical potential is chosen to 
yield 16% doping (/i = — i). 

The local impurity is introduced into Hamiltonian 
Eq. (2) by modifying the electron energy on a particular 
site [Til jj. The corresponding addition to the Hamil- 
tonian is 



#imp = Mmp(«0T + n oi) + Si m p(n 0T - n i). 



(3) 



The first term is the potential part of the impurity energy 
that couples to the total electronic density on site 0, and 
the second term describes the magnetic interaction of the 
impurity spin and the electronic spin density on the same 
site. We assume that the impurity spin is large and can 
be treated classically, as if it were a local magnetic field. 
The goal is to determine Vi mp and Si mp so as to match 
both the location of the impurity states within the gap 
and the spatial distribution of their intensity. 

First let us analyze the position of the impurity energy 
level as a function of the impurity potential. We solve the 
MF equations self-consistently on a square lattice with 
periodic boundary conditions. The results are presented 
in Figure ^|. All impurities placed in a superconductor 
generate quasiparticle weight both on positive and sym- 
metric negative biases. The sign of the energy level is 
defined based on where the majority of the quasiparti- 
cle weight resides. Attractive impurities produce energy 
levels lying below the Fermi surface, while repulsive ones 
generate states with positive energy. In the limit of a 
very strong impurity potential both attractive and re- 
pulsive impurities generate identical states with a small 
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residual energy related to the amount of the particle-hole 
symmetry breaking, caused by the specifics of the band 
structure and doping p2[. 
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FIG. 2. Impurity energy level position as a function of the 
potential impurity strength, Vj mp . The blue (top) line corre- 
sponds to repulsive impurities, which produce positive energy 
levels (eimp > 0); the red line is for attractive impurities, 
which generate negative energy states (ei mp < 0). For strong 
impurity potential both energies converge to the same small 
value determined by the amount of particle-hole symmetry 
breaking. 



The analysis of the energy level positions implies that 
the Zn impurity can be associated with a strong at- 
tractive potential, Vz n = —Hi = —4.4 eV and ezn = 
— 0.005< = —2 mcV. Ignoring for now the level splitting, 
the Ni case can be associated with a relatively weak re- 
pulsive impurity, VNi = t and en; = 0.0443£ = 18 meV. 
These impurity strengths are in agreement with the gen- 
eral band structure arguments. Zn 2+ ion has 10 electrons 
that completely fill d orbitals. Hence, the d x 2_ y 2 orbital 
of Zn, relevant for interaction with Cu-0 plane orbitals, 
is deep below the Fermi surface. On the other hand, 
Ni 2+ has 8 electrons in the d shell, with the d x 2_ y 2 being 
unoccupied, but with a small energy, given by the level 
splitting within the d shell. 

The spatial distribution of the spectral intensity corre- 
sponding to the Zn impurity is shown in Figure |[ The 
top two plots show the intensity Af™° as it would be seen 
if STM tip were directly imaging DOS in the Cu-0 layer. 
The intensity is related to impurity state wavefunction 'F 
on site (i,j), 
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The intensity on the impurity site is suppressed due to 
the strong impurity potential. The bottom two plots cor- 
respond to imaging through the top Bi-0 layer. They are 
obtained by applying a filtering function F[^] = \Mij\ 2 



to the impurity state wavefunction '5. The effect of the 
filtering function is to produce the intensity 
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Indeed, the intensity on the Bi-0 layer is maximized on 
the impurity site due to the interference of the contribu- 
tions from the impurity's nearest neighbor sites. 
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FIG. 3. Top two plots show the real-space spectral inten- 
sity of the calculated Zn impurity state in the Cu-0 plane. 
While most of the intensity is concentrated at the negative 
bias, -2 meV, as it is seen in the experiments Q, the spatial 
shape of the state does not agree with the experimental re- 
sults. The bottom two plots show the same impurity state, 
but as seen through the Bi-0 layer. The figure is obtained 
by applying the filter function of Eq. (Q). Both the energy 
of the state and the spatial intensity distribution agree with 
experiment jt|. 

The structure of the Ni-induced state is more compli- 
cated than the Zn case. It appears on both positive and 
negative biases. In addition, there is a peak splitting on 
each bias. Unlike Zn 2+ , Ni 2+ impurity is magnetic, with 
spin 1. To simulate the effect of spin we include a non- 
zero magnetic part of the impurity potential, 5i mp 0, in 
the Hamiltonian Eq. (||) . The spin component introduces 
level splitting between "up" and "down" spin states. Fig- 
ure ||a shows the amplitudes of the spin-split states on 
the impurity sites and its neighbors for Vm = t and 
Sm = 0.4t. The spin-split energy levels are ej_ = 0.052< 
and £| = 0.037i. The total spin-up and spin-down inten- 
sities for each bias are shown in Figure gb. Both fi gures 
(a) and (b) correspond to the image affected by the fork 
filter. The general shape of the states agrees well with 
the experimental data j^j . The total weight on the nega- 
tive bias is about factor of 3 smaller than the weight on 
the positive bias. 
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FIG. 4. Simulation of Ni as an impurity with weak mixed 
potential, Vm = £ and SW = 0.4t. (a) Spectral intensity on 
the impurity site, its nearest neighbors (nn), and next nearest 
neighbor (nnn). Red line is spin up, and blue is spin down, 
(b) Intensity map for the combined spin up and spin down 
states. 



Some aspects of the Ni-impurity state cannot be ad- 
dressed within the simple framework of our theory. We 
assumed that the spin of the impurity is static and plays 
the role of a local magnetic field. Bulk measurements in- 
dicate that Ni dopants do not introduce additional mag- 
netic moments. This may be either due to the partial 
Kondo screening of S = 1 Ni spin or due to the local 
attraction of additional hole ("self-doping") by Ni 2+ to 
form an effective Ni 3+ complex with spin 1/2. No mag- 
netic signature, combined with the weak potential char- 
acter of Ni may explain the relatively weak pair-braking 
effect of Ni doping observed experimentally 

In the case of Zn, we focused on the large local po- 
tential scattering off the impurity as a main effect to 
account for the STM intensity distribution. It is impor- 
tant to mention that Zn substitution also induces local 
magnetic moment, which is observed in magnetic sus- 
ceptibility and NMR measurements jl3| . The quenching 
of magnetic moment around Zn can lead to the Kondo 
resonance near the Fermi surface, as has been proposed 
by Polkovnikov et al. JL4J. Below the Kondo tempera- 



ture, Zn should behave as a unitary scatterer, virtually 
indistinguishable from either local or extended |l4|,[l5|] po- 
tential scatterer. We argue however, that the tunneling 
fork filter should apply regardless of the details of the 
scattering mechanism. 

The tunneling fork mechanism has consequences be- 
yond the settings of the original experiments ■ Here 
we point out two of them: (1) the direct tunneling into 
Cu-0 plane, and (2) the structure of the Cu vacancy 
states. One could assume that the tunneling into the 
exposed top Cu-0 plane should be free of the tunnel- 
ing fork. However, upon closer inspection it is clear that 
even in this case, the direct tunneling from tip into the 
planar 3d x 2_ y 2 orbitals is exponentially weak compared 
to the indirect tunneling through the s-wavc like orbitals 
of Cu, Zn, or Ni extending out of the plane. Hence, the 
tunneling fork should apply also in this case, resulting 
in the same spatial form of the filtered impurity states 
as obtained above. In the case of Cu vacancy in the 
Cu-0 plane, under assumption that there are no s-wave- 
like orbitals centered at the vacancy site, the tunneling 
fork mechanism implies that resonance state observed by 
STM should be the same as for Zn, but with no spectral 
intensity in the center of the pattern. 

In conclusion, we have demonstrated that recently ob- 
served Zn and Ni impurity states in BSCCO can 
be explained by simple model of strong potential im- 
purity in the case of Zn and mixed (potential + spin) 
impurity in the case of Ni, interacting with a d-wave su- 
perconducting condensate. The crucial aspect that we 
have included in the present treatment is the effect of 
the quantum-mechanical tunneling between the STM tip 
and the Cu-0 planes. The same non-local filter effect 
should be operable in the case of the "direct" tunneling 
into Cu-0 planes and should have consequences for other 
types of defects, e.g. Cu vacancies. 
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